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Spin relaxation in a site-disordered Ising model within master equa- 
tion approach is studied. The q, w-dependent susceptibility of the 
OO ■ model is calculated and investigated. Effects described by the two- 

site cluster approximation and lost by the mean field approximation 
are discussed. Comparison of obtained results with dielectric measure- 
ments in Cs{Hi^xDx)2P04 is presented. 

> 

^ ' 1. Introduction 

I ! 

\^ . Being the model of critical phenomena in many different objects an Ising 

' model compels an attention during many years. This paper investigate the 

I kinetics of the model (within master equation [ ||] approach) having sev- 

■ eral goals. We consider the disordered (site chaos) Ising model in order 

to describe the dynamic properties of partially deuterrated ferroelectrics 
with hydrogen bonds. Besides the theory seems applicable to many other 
systems: alloys of magnets, diluted magnetic or ferroelectric systems and 
others. Within two-site cluster approximation (TCA) we calculate the dy- 
namic stucture factor to supplement our previous study [ ^ of correlation 
functions of the model. 



2. The Hamiltonian. The approach 



^ ' We consider the spin system with site chaos on the Bravais lattice which is 



described by the Ising Hamiltonian 



H=-Y.KiSi-\Yl KijSiSj ■ (2.1) 



H describes the set of Ising spins with the pair exchange interaction Kij in 
the site-dependent field Kj. The set of variables {Si} {Si = ±1; i = 1 ■ ■ ■ N, 
N is a number of the sites) represents a state of spin subsystem, and the 
set {Xia} describes the sort configuration {Xia = 1, if the site i is occupied 
by spin of the sort a, otherwise Xia = 0; a = 1 ■ ■ ■ Q, Q is a number of the 
sorts). The quantities Ki, Kij depend on sort configuration 



l^iaXia', Kij — ^ ^ KiQi^jpXicfXjp (2-2) 



^ ^ ■ ^ l/t^ fcLX) L J / J 
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SO the true parameters of the Hamiltonian are Kia, Kiajfs- Here we 

suppose the pair interactions to be short-range 

f 1 if 7 G TT' 

Kia,jf3 = KafiTTij-, K^p > 0; TTij = <^ q] otherwise (^'^^ 

where vTj denotes the set of the nearest neighbours of the site i (the first 
coordination sphere). 

From the very beginning we shah consider the quenched system (frozen 
sort configuration) and investigate its dynamics on the basis of Glauber's 
master equation [ ||] for the distribution function p{{S} ,t): 

j^Pi{S},t) = -J2W,{---Sr--)p{---Sr--,t) 

i 

+ Y.Wi{----Si---)p{----Si---,t) (2.4) 

i 

Following ref. [ ^ the probability per unit time that the ith spin fiips from 
the value Si to —Si is assumed to be 

Wi{---Si---) = -3-(l-5,tanh/?e,) (2.5) 

iTi 

where 

ei = Ki + Y^KijSj (2.6) 

i 

is a local effective field acting on the ith. spin. In this approach the dis- 
sipative properties of the spin subsystem is caused by its interaction with 
another (e.g. phonon) subsystems which play the role of a heat bath. The 
quantity Ti has the dimension of time and defines the time scale of relax- 
ational phenomena. We assume the quantity to be locally sort-dependent: 

n = Y,rlX,^ (2.7) 

a 

Relations ( |2.4D and (^), ( |2.6| ) lead to the following kinetic equations for 
the average value of spin: 

Di^t{Si)H,t = {i^^^Pei)H,t (2.8) 

where 

{{■■■))H,t = Spsp{{S],t){---); D,^t = l + Tij^ (2.9) 

The system is disordered, therefore in order to calculate experimentally 
measured quantities one have to perform an averaging over the sort config- 
urations. For example, the polariza tion of a partially deuterated H-bonded 
ferroelectric with the Hamiltonian (2.1) is defined as 

^(*) = ^Y.Y.^'o.{{Sio^)H,t)x; Sia = SiXia (2.10) 
i a 

where y is a volume of the system, //q the dipole electric moment associ- 
ated with the quasispin of sort a (proton or deuteron), {■ ■ ■)x denotes the 
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averaging over sort configurations with some density matrix p{{Xia}) that 
depends on conditions of system's freezing. The two-site cluster approxi- 
mation, which is used below, is sensitive only to the following moments of 
the distribution p{{Xia}) 

{Xia)x = Ca, {XiaXjp)^ = Wa/S (j G T^i) (2.11) 

where Ca has the meaning of the sort a spins' concentration. 
3. Formulation of two-site cluster approximation 



The equation (2^) is quite intractable. One can expand the {tanh. Pei)H,t in 



the right-hand side of (2.8) and find that this equation couples correlation 
functions (CFs) of the type {Si^ •■■Si^)H,t where the site indices belong 
to the first coordination sphere of the site i ( ii ■ ■ - in G T^i)- The equation 
for {Si^ ■ ■ ■ Si^)H,t derived from ( [2.4]) and ( |2.5| ) involves higher CFs, so the 
exact treatment of the Glauber's equations has succeded only in the case of 
one-dimensional system in zero field (ideal (one-sort) chain [|l|, ^, ^, chain 
with single impurity [|5|). In other cases interpolative approximations are 
used which express higher CFs via lower, e.g. (Hj Sj)H,t ~ Y\j{Sj)H,t- 

We suggest here the closure to the equation ( |2.8| ) in the spirit of a cluster 
approach. In a one-site approximation we replace contributions of all spins 
with effective fields 

where r'fia has the meaning of effective field acting on the spin i of sort 
a from the nearest neighbour on the site r. In a two-site approximation a 
contribution of one of the nearest neighbours is taken into account explicitly 

Ei ef^ = Ki^t + + ^ij^j = j^i't + KijSj, {j G TTj) . (3.2) 

r / 2 

In the case of ideal (one-sort) system no more assumptions are needed to 
obtain the closed set of equations for the average value of spin {Si)H,t 

Di,t{Si)H,t = {tanh pef})H,t . (3.3) 

For the many-sort system it is necessary to point out a dependence of the 
fields r'-Pi,t on the sort configuration. We shall believe the fields r'fu to be 
different in the one-site approximations and the two-site one until averaging 
over sort configurations is carried out. For the sort averaging procedure both 
"one-site" and "two-site" fields will be assumed to be locally sort-dependent 

a 

where r'fia,t is a c- number (it does not contain spin or sort operators). 
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After simple transformation the equation ( p. 3D can be written in the 
form 

Di^t{Si)H,t = tanh/5Kj,i 

Di^t{Si)H,t = Lij^t + Pij,t{Sj)H,t (3.5) 
Dj,t{Sj)H,t = Lji^t + Pji,t{Si)H,u 3 £ TTj) 



where 



Lij^t = \ [tanh^( jKjf + Kij) + tanh/3( jKjt - K-ij)] 

Pij,t = \ [tanh/3( jKjf + Kij) - iaxih. {3 { jRu - Kij)] (3.6) 



The expressions (3.6) consist a system of nonhnear differential equations. 
To linearize this expressions we consider only linear responce of the system 
on a small time-dependent changing of the external field. For this purpose 
we write the dynamic quantities in the following form: 



l^i,t — f^i ~\~ 5f^i,t] f^i,t — f^i ~\~ ^Ki^t') 

j^i,t = jKi + 6 jRi^t; {Si)H,t = "T'f ^ + '^"T'l.V; 

-^ij,t — Lij + Lij 5 jRi^t'i Pij,t — Pij + Pij ^ j^i,U 

tanh/3Kj_t = tanhKj + (1 — tanh'^ k,) 5Ri^t (3.7) 

Here Kj, Ri, jRi, m[^\ Lij, Pij are the static (independent of time) com- 
ponents of corresponding quantities and 6Ri^t, ^ j^i,t, ^fn'it ^ ^'ij ^ j^i,t, 
Pij 6 jRi^t, represent linear in 5Ki^t deviations of these quantities from their 
static values. 



Lij = - (tanh/?( + Kij) + tanh/?( - Kij)) 
Pij = t; {tanhP{ jRi + Kij) - tanhp{ jRi - Kij)) 



^ij = ^ (cosh ^ PijKi + Kij) + cosh ^ PijKi - Kij] 

Pij = ^ (cosh-2 p{jR, + Kij) - cosh-2 Pi^Ri - K,j)) . (3.8) 

Due to ( t3.7D two independent system of equations is obtained from ( ^.61) . 
The static part: 

mf^=tanh/3Ki (3.9) 
mf ^ = Lij + Pijmp; mf^ = Lji + Pjivaf^ 
can be reduced to the form 

= tanh/?K, = y . (3.10) 

The same equation is also obtained on the basis of cluster expansion of the 
free energy [ |6| and can be written in the form 



(1) 



{S,),^ = {S,)p^^ (3.11) 
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where the one-site pi and intracluster pij density matrices are defined as 
follows 

_ exp(-/3gO TT _ - Q 

Pi — T< 7 ~ 

bps^ exp{-pHi) 

Pij — 7; 7 TTTr 7' -^ij ~ j^i^i ~l" i^jSj + KijSiSj (3.12) 

Sps^Sj^y^v{-PHij) 

Multiplying (|3lo| ) by and averaging obtained relation over the sort 
configurations we get the equation for an average value of sort a spins 

"^ia = tanh Ria 

^ smh.(3{jRia + iRjfi) + aapsmh.[3{jKia - iRjp) ,„ 

= / — —. — z = ^ — 777 — Z Z T (o.ioj 

^ cosh (3 [jKia + iKjf3) + UafS cosh fJ[jKia - iKjp) 

On assuming uniform field — > Kq (UF: rfia — ^ j^ia — ^ ~l~ — 
l)(pa = R'a, i^ia K,a+z(pa = Ra, z\s a first Coordination number) we obtain 
the known in corresponding static theories (TCA [|2|, first order of cluster 
variation method [ ^) relation: ( 3.13| ) becomes the set of 20, equations for 

the same number of quantities - order parameters rrSa^ and $7 effective 
fields (pa ■ 



The dynamic linearized part of ( |3.6| ) (here nonlinear on 5Ki^t terms 
L'- 5rni^ SiRj^t, L'j^ 6m^^^ ^ jRi,t should be neglected) after time Fourier 
transformation and averaging over the sort configurations yields the Orn- 
stein-Zernike-type equation for the dynamic structure factor i^^lJ id^'^) ~ 



(1) ^'^''^^ 
^SikI (uj) '-'^ model. Spatial Fourier transformation presents its solution 
in the form 



+ (^g - .,) f,^(2°) U^^'A'' <^(2°) - ) (3.14) 



where the matrix elements are 



(2)\ UF ^ l-tanh^/3Ko 



Jaf3 "^4- DaD^ - Pa-yPya 
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UF R/3aPaf3 



(m(2)(g,a;)) = E ^''■^'^^"'^^^"^S,73(^) (3-15) 



'a/3 
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and 



a/3 



Rap 



- (tanh/3(K^ + K^p) - tanh/3(K'„ - K^p)) (3.16) 



2a. 



aP 



2aap{^ + cosh 2/3k^ cosh2/?K^)+ 



(1 + alp) {cosh 2 l3R'a + cosh2/3K^ 



(1 + alp + 2aap cosh 2/3k^ 



cosh + kJj) + a^/j cosh 



Factor vr^- has the form 



(3.17) 



Annealed model 



In this case protons (deuterons) can go from one site (hydrogen bond) to 
another. Therefore the sort configuration changes in time. This motion 
should be described by Kawasaki master equation But we assume that 
the sort configuration changes slowly and for the frequencies of magnetic (or 
ferroelectric) dispersion quasispins see sort configuration at rest. Therefore 
the formula ( ^.14 ) remains valid. If the temperature changes quite slowly, 
Waip is a temperature-dependent quantity which can be found by the static 
theory [ |]. 



5. Discussion. Numerical results. 



The expression (|3.14D is our target formula. Let us consider its partial cases. 
All expressions will be given for UF case. 

1) Noninteracting system {Kap = Vq, /3) is a trivial case. 



m 



= Ca tanh 



''^apil^^) = ^^13 1 



Ca(l — tanh^ P^a] 







(5.1) 
(5.2) 



All the spins are independent, so the structure factor are independent of 
the wavelength. The relaxation times are equal to r^. 

2) One-sort (ideal) system. This case can be obtained in the limit ci = 

1. The single order parameter m^^^ and the effective field (pi satisfy the 
following equation 



m«/ci = cx«=tanh/3^,- ^^"^'^^'i 



where 



cosh 2/3 k'^ + ail 



Ki = Ki + zipi] k'i = ki + zVi) z' = z — 1; a = e 



-2K 



(5.3) 
(5.4) 
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f'2'1 

The dynamic structure factor m\j^ {q, u) takes the form 

^11 = ni' {q, uj) = ^ ^ .^^^^ , (5.5) 



Figure 1: Pure system on plane square lattice (z = 4). Static pair CF 

(2) -» 

m\^ {q,0), inverse relaxation time l/r(q), real and imaginary parts of dy- 

namic pair CF m\^{q,u:i) at different wavevectors: VTq* equals 4 (g = 0), 2, 
0, —2, —4 (the last corresponds to the vertices of the first Brillouin zone). 



where 
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F 



l + a^ + 2a cosh 2(5R' ^ _ / (i) 



and Gil {q, 0) coincides with a spatial Fourier-transform of the static corre- 
lation function {SiSj)'jj, which was founded with cluster expansion of the 
free energy within TCA in the ref. [ ^. Result ( p. 61 ) is accurate in the 
one-dimensional case. Relaxation time 

of, P Y^T>T, of, _ tanh/?K 



T{q) = tH 1 - vr„- -—^ '-^^ r'^ 1 - 7r„- ^ , (5.7) 

^' H ^l + z'P2; 1^ ''l + z'tanh^/?^; ^ ' 

is a decreasing function of q. Note, that known exact result (zero external 
field, ideal one-dimensional system) is a partial case of the formula ( |5.5D . 
Figure 1 demonstrates the temperature and g-dependence of the static 

structure factor m^/(^, 0) and the relaxation time r(g) of the system on the 
plane square lattice (z = 4). Five lines represent five wavevectors: tx^ = 4, 
2, 0, —2, —4. In this figure (as well as in the formula (^.T])) one can see an 
interesting symmetry of inverse relaxation time: tau\/tau{q) line at tt^ is 
mirror reflection of the line at —tt^. The figure 1 shows also a behaviour of 
real and imaginary parts of the dynamic structure factor (at the frequency 
oj = O.l/r?). 

3. Diluted system (interacting spins (sort 1) are diluted by the noninter- 
acting (nonmagnetic, nonferroelectric) impurities: Kia = for all a ^ 1) is 
interesting due to percolation phenomena in such system: below some con- 
centration of sort 1 spins there is no infinite cluster of interacting spins and 
ordered phase does not appear. The mean field approximation (MFA) fails 
to describe percolation phenomena [ |^. At zero external field within MFA 

and pair CF per interacting spin depend on temperature, concen- 
tration and interaction strength only via T/T^^^ [T^^^ = zciKu/ki^): 
ai^) = cjf)(r/r5^^), 4? = a^^^{q,u;,T/T^^^^). Therefore MFA predicts 
that (7^^ reaches its saturation value {cr^^ = 1 at T = 0) at any concen- 

(2) 

tration of interacting spins, and m\i {q,Lo) shows a single relaxation time 
behavior. In TCA = ^ , 

(I) , sinh2/?K'i , 

m\ = ci tanhpKi = wu — - — — ; 1- (ci — wn) tanhp/ti (5.8) 

cosh2pK;^ + ail 

and this takes place for the case wia = only, otherwise at zero temperature 
a\^^ < 1 (due to the presence of finite size clusters of interacting spins) and 
becomes zero if ci > z'wu (below percolation point). In the case of complete 
chaos {wa/s = CaCp) TCA predicts percolation concentration Cp = [z — 1)^^. 
We note that the Effective Field Approximation of Kaneyoshi et.al. [ ^ 
(magnetization and static susceptibility) also reproduces this property of 
the diluted system. 

Dynamic pair CF contains three terms 
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where 



T>Tc C1W12 
° ZCiFf ) - Z'W12R12 ^1 + ^'^11 
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Figure 2: Diluted system (ci = 0.6, wu = cf) on plane square lattice 
(percolation concctration Cp = 1/3). Amplitudes A^, A_, Aq, inverse relax- 

ation times l/r_|_(g), l/r_(g), and dynamic pair CF m\i{q,u)) at different 
q: 7r,-=4,2,0,-2,-4. 
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Figure 3: Diluted system. Temperature and frequency dependence of dy- 

(2) 

namic pair CF m\i {q, u) at q = 0. 



and 

= wiiRii + W12R12 



2{Ri) (^ciFf ) - z'{Ri)) 

V{Tr^) = (^ciFi'^^wnRmr^y 

+ A{Ri) (zciFf ^ - z'{Ri)) W12R12 (^ciFf ^ - /i2i2)(5.11) 
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Figure 4: Spontaneous polarization P, inverse permittivity 1/e and inverse 
relaxation time 1/r of the ferroelectric C s{Hi-xDx)2P04 for different de- 
grees of deuteration. Squares correspond to experimental data of [ 12|, bold 



lines on the lower picture connect experimental points from the ref. [ 13|. 
Thick lines ~ theoretical results for annealed system with parameters given 
in table 1 and concentrations a; = 0, 0.18, 0.54, 0.88, 1. 
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Figure 5: Comparison of experiments and quenched model with assumption 
Wag = CaCg (complete chaos) (see caption to the previous figure). Thick 
lines - theoretical results for quenched system with parameters given in 
table 1 and concentrations x = 0, 0.12, 0.47, 0.83, 1. 
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quenched annealed 


Kii/kB, K 




610 


Ki2/kB, K 




500 500 


K22/kB, K 




390 


Ju/kB, K 




2.85 


Ju/kB, -ft- 




3.2 1.6 


Jon/kn K 
'J 221 i^Bi 




l.Ob 


{Vii + V22 - 


-2Vi2)/kB, K 


- -1000 


lafs/kB, K 







/ii, e.s.u. X 


cm 


2.18x10"^^ 


^2, e.s.u. X 


cm 


1.87x10"^^ 






3x10-14 






5.5x10-1^ 



Table 1: The model parameters describing experimental data on the 
pseudo-one-dimensional ferroelectric Cs{Hi-xDx)2P04: short-range cou- 
pling Ka/3; long-range interaction Jq,^; effective dipole electric moments //q 
associated with proton or deuteron on hydrogen bond; parameters of mas- 
ter equation - relaxation times of noninteracting protons or deuterons. 
Vai3, Iaf3 are parameters describing nonexchange interactions between pseu- 
dospins [2]. Deuterons correspond to the pseudospins of sort 1. 



Ff) = c,(l-(a«)^) 
R12 = 1 — tanh^ Pk'i 

1 + a cosh 2Pr'i 



(cosh2/3K; + a){l + a'^ + 2acosh2(3R[] 



i [tanh P{k[ + K)- tanh I3{r[ - K)] (5.12) 



The first term of ( ^.91) has the (^.2|)-like behavior (relaxation time equals 
r}*), so we believe it to be a contribution of isolated spins (i.e. the spins, 
which all neighbours are noninteracting). The last two terms of ( |5.9| ) have 
the following symmetry 

B±{ti^) = B^{-7T^y, A+{-z) = A^{z) = 

D±{7rg') = -D^{-TT^); D+{7i^) > 0; D_(7r,-) < (5.13) 

The relaxation time t_ is always less then rf* while r+ is always greater then 

Ti and slows down at the critical point (r+(0) — >' 00, A^{'Kfj) — 00). 
Both r_ and increase with the approaching to the center of the first 
Brillouin zone (g ^ 0). The amplitude A- is equal to in the centre of the 
first Brillouin zone {iTg- = —z) while Aj^ vanishes at = —z (in the vertices 
of the first Brillouin zone). 

In the case wn = ci {wia = Va / 1, all interacting spins belong to 
one infinite cluster), Aq is equal to zero. A- equals zero at -K^q) > and 

A^ is equal to zero at 7r(g) < 0; cj[^^ and pair CF are given by the formulae 

(HI, (El- 
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Figure 2 shows the pair CF (|5.9D of the diluted system on the plane 
square lattice (z = 4). In the figure 3 one can observe that presence of many 
relaxation times in dynamic pair CF of the diluted system can obviously 
be seen at small temperatures and frequencies. Abscissas of maxima of 

— Im m^ii{0,uj) approximately equal inverse relaxation times: ujq ~ l/''"?) 

~* (2) ~* 

ujj^ « 1/t_|_(0). Specific temperature dependence of m]^/(0,u;) is connected 
with nonmonotonouos behavior of T-|-(0). 



6. Experiments 

One can see that TCA reflects essential properties of the diluted system. It 
seems to us that quality of the approximation is even higher for the system 
without noninteracting impurities {K^p 7^ 0) because such "effective field" 
theories are more appropriate in the case of small fluctuations in a system 
considered. Obtained within TCA characteristics can be used in order to de- 
scribe a variety of experimental data on partially deuterrated ferroelectrics 
with hydrogen bonds and alloys of magnets. For example, good fit of TCA 
results for the ideal (one-sort) system to experimental data on the ferroelec- 
tric dispersion in pseudo-one-dimensional {z = 2) ferroelectrics CSH2PO4 
and its deuterated analog CSD2PO4 has been reached in the ref. [ |ll[|. 
The many-sort theory developed here is applicable to partially deuterated 
compound - Cs{Hi-xDx)2P0/^. We obtain a good agreement of TCA re- 



sults on annealed model and results of dielectric measurements [12, 13] 
(figure 4). Fit of the quenched model (figure 5) to experimental data is not 
so good. The model parameters are shown in the table 1. More information 
about an actual type of the system can be given by measurement in the 
low-temperature region, where significant differences between behaviours of 
the annealed and quenched models exist, and by the wavevector-dependent 
experiments (e.g., neutron scattering). 
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Appendix 

Here we consider a system with the Hamiltonian 

— H/ksT = - ^ KiajpSiaSjp + - ^ Jia,jl3SiaSjp + 'Y^^iaSia , (Al) 

where Kiajjs, Jiajp are pair interaction parameters, Tia external field acting 
on the spin of sort a on the site i. The interaction Ki^jj^ is same as in (2.1)- 
(2.3), but Jia,ji3 is assumed to be of long-range type, i.e. not restricted to 
the nearest neighbour sites. The Hamiltonian (Al) can be written as 



n- H — Jiajp'^ A m-p (A2j 

ia,jf3 

MFAtt _ TT ,^bT ^ La)Lil) 



ia,jf3 

r 

2 ~ 

ia,jf3 



where H is the Hamiltonian ( p.l| ) in which the field k is taken in the form 
MFA assumes 

^mSi^ = {{S,^)n,t). (A4) 
and ignores last term of (|A2|): 

-ly , MFA u _ u , >^ Ln)Ln) /acn 

H — > H - H + 2^ Jio,jp m^^ m-f^ . (A5) 

ia,jl3 

One can see that formalism of the section 3 is applicable to the system with 
the Hamiltonian (|A5|). The equation of motion has the same form 

Di{Si)H,t = (tanh/3ej)7^^j; = ^KijSj + m , (A6) 



DYNAMICS OF DISORDERED ISING MODEL 



17 



but now the field includes molecular field (see (A3)) and thus couples the 

(A7) 



selfconsistency parameter ^m[l^ . Then one can find that 



L (1) (1) 



where m-^^ is defined by the equation (|3.13D and 



(SijS^^ + J2Jiy,k5''^klpi^)] • (A8) 

\ kS ) 



(A8) leads to the following expression for the structure factor 
(^m(2)(g-,^))"' = (rn(2)(^-,a.))"' - J{q) , 



(A9) 



where i7i^'^\q,iv) is defined by ( 3.14 ) and (j{q)) is a Fourier-transform of 

V / a/3 

Jia,j/3 ■ 
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